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Abstract 

We reanalyze high energy QCD scattering regimes from scat- 
tering in cut-off AdS via gravity-gauge dualities (a la Polchinski- 
Strassler). We look at 't Hooft scattering, Regge behaviour and 
black hole creation in AdS. Black hole creation in the gravity dual 
is analyzed via gravitational Shockwave collisions. We prove the 
saturation of the QCD Froissart unitarity bound, corresponding to 
the creation of black holes of AdS size, as suggested by Giddings. 



1 



1 Introduction 



Gravity-gauge dualities have been an important tool in getting information about nonper- 
turbative Yang-Mills theories, since the original work of Maldacena (see also [2]). But the 
contact with real high energy experiments has been lacking, partly because of the absence 
of a gravity dual of QCD. However, in the work of Polchinski and Strassler (and the 
later [3] dealing with deep inelastic scattering) it was shown that one can derive a lot of 
information from just a simple gravity dual model of AdS cut-off in the IR (where the mass 
gap is related to a modification of the geometry) and maybe in the UV, in other words the 
2-brane Randall- Sundrum model [S]. 

High energy QCD scattering of colourless objects (e.g. glueballs) is thus related to scat- 
tering in the cut-off AdS via convoluting the scattering amplitude with AdS wavefunctions. 
This simple model still allows one to get a lot of information about QCD at high energy, 
when the IR modifications due to the mass gap are not so important. 

Giddings [6] took this proposal further and analyzed the scattering inside AdS in the 
extreme inelastic case, when black holes are being formed. He analyzed the black hole 
formation using a simple model, devised in [HE] (see also earlier work in [5]), the cross section 
for black hole formation being equated with the geometric area of the black hole horizon of 
mass equal to the center of mass energy (er = ixr 2 H , r H = Th(\/s))- At moderate energies, 
the black holes being formed can be taken to be in flat space, but at higher energies, the 
size of AdS becomes important. Giddings argued, by calculating the Newton-like potential 
(linearized Einstein gravity), that the horizon of a black hole formed on the RS IR brane 
will be such as to give a saturation of the Froissart unitarity bound for the cross section: 

e -Afir i i 

h 0;lin ~ G 4 y/s => r H ~ — ln(G 4 y/sM^ => a ~ irr 2 H ~ (xrMv^)) 2 (1.1) 

r Mi Mi 

But there are a lot of uncertainties about this calculation. In a previous paper |lUj . 
we addressed some of these uncertainties, by looking at black hole formation via shockwave 
collisions (a technique first proposed in and estimated the maximum impact parameter 
bmax for which a black hole is being formed, and calculated the minimum mass of the formed 
black hole. The calculation was done in flat D > 4, as well as in the background of the one 
brane RS model ("alternative to compactification" [12] )> as was apropriate for the various 
low Mpi scenarios ^3 El E) in which one could detect black holes at accelerators. We have 
set up the general formalism in a way to be used for the AdS and 2-brane RS calculations 
needed for the QCD dual scattering. We also analyzed the effect of string (a') corrections 
to the black hole creation. Quantum corrections were also analyzed in a different way in 
|15| ITo] and 't Hooft scattering inside AdS was analyzed using different methods in [TJj. The 
advantage in this semiclassical formalism is that now a = Kb 2 max is rigorous. 

In this paper we will explore the consequences of the calculations in ^U] fo r the high 
energy QCD scattering, and we will revisit some of the previous results to see if we can gain 
more insight. We will try to use consistently the Polchinski-Strassler set-up, in particular 
finding how to turn the classical scattering with black hole formation, happening at a certain 
point in the gravity dual, into an integration over the gravity dual. We will use a simple 
black disk model to turn the classical scattering into an imaginary elastic amplitude for 
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which we can use the P-S formalism. The most important piece of information learned in 
this paper will be that we are able to justify in a more rigorous way the appearence of 
the Froissart bound. Given the Polchinski- Strassler set-up, we will calculate the maximum 
impact parameter being formed in the scattering of two Shockwaves, and thus get a = ^b 2 max . 
The importance of this formalism is that the Shockwaves are exact solutions (not linearized 
ones), giving an advantage over the horizon calculation in jB]. 

We should note here that we are not attempting to define a general amplitude (nor a 
S matrix) for a general scattering inside AdS. We will see in the next two subsections that 
for our purposes it is enough to assume the existence of a nonperturbative amplitude (and 
corresponding cross section) when the external legs are situated in a flat 4d slice of AdS 
space. We will define a simple black disk model for this amplitude that mimics the (as yet 
unknown and undefined) nonperturbative amplitude by giving the same cross section. For 
external legs living in a flat 4d slice, one can certainly define amplitudes and cross sections 
in the usual manner. 

The paper is organized as follows. In section 2 we review the formalism of high energy 
scattering in QCD from scattering in cut-off AdS. In section 2.1 we show how to take a 
classical scattering in the gravity dual, characterized by a b max (s), and turn it into a quantum 
amplitude, which we can relate to QCD. In section 3 we describe the general formalism for 
calculating b max (s) for black hole formation in the scattering of two Aichelburg-Sexl (A-S) 
waves in a background of AdS type. In section 4 we apply this formalism for the case 
of scattering inside AdS and on the IR brane in the 2-brane RS model. Section 5 is the 
most important, putting together all the pieces of information, and analyzing the various 
QCD energy regimes, and discussing the saturation of the Froissart bound, as well as string 
corrections in the gravity dual. In section 6 we conclude and in the Appendix we show the 
details of the calculation of the trapped surfaces being created when two A-S Shockwaves 
scatter inside AdS and on the IR brane. 



2 High energy QCD from AdS 

Polchinski and Strassler jHj have found a simple model for relating high energy QCD scat- 
tering with scattering inside AdS. 

For a conformal field theory, the corresponding near horizon (r — > 0) metric for the brane 
configuration is 

ds 2 = ^da? + %dr 2 + R 2 ds 2 x (2.1) 
Rr r z 

The global momentum = —id^ (momentum for gauge theory scattering, for instance) is 
thus related to the local inertial momentum (of a local inertial observer in AdS) by 

% = jp» (p^tT = PvPvgT) (2.2) 

Thus high energy is large r and low energy is small r. And then a nonconformal gauge 
theory like QCD will just be modified at small r (low energy), at high energy remaining 
conformal. 
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Since the low energy cutoff for the conformality of the gauge theory should be of the 
order of the mass of the lightest glueball state A, the cutoff on the AdS geometry is 
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R 2 A (2.3) 



Thus in theories with mass gap like QCD the warp factor becomes bounded, and the simplest 
effective description that we are going to use throughout the paper is just to cut-off the 
integration over r at r m j„. We can also cut off the theory at high energy, and with UV and 
IR cut-offs we have the 2-brane Randall-Sundrum model [Sj. 

Corresponding to the string tension in AdS there is also a gauge theory string tension 
a' = {gN)- l / 2 A7 2 and 

/ 7~ nr. r min / 7 _ , nr. In A\ 

yap = yap > yap < yap (2-4) 

r 

A glueball correponds in AdS to a state with wavefunction 

e ipx i>(r,n) (2.5) 

(plane wave in 4d and some wavefunction for r and X). Assuming that scattering of gauge 
invariant states (e.g. glueballs) within Yang-Mills is equated by AdS-CFT with a scattering 
inside AdS of the above states; moreover assuming that the states scatter locally according 
to the flat space amplitude, we get 



<Mjp) = J drdrQy/gAatri-gip) IJ^ 



(2.6) 



Then integrating over r corresponds to integrating over scattering energies in the local frame 
(p). In this picture a' is a constant, as in flat space string theory, but momenta (p) "run" 
as r is varied, as can be seen from (|2.4|) . It is also clear from (|2.4J) that one can consider p 
to be fixed and a' to run, but our interpretation makes it clear that we are integrating over 
string theory momenta. One should also note that although originally in the formula ()2.6j) 
A string (p) was meant as a worldsheet string amplitude (with vertex operator insertions), that 
in flat space gives the perturbative spacetime amplitude for scattering particle states (e.g. 
gravitons), we are now extending its meaning. We will first consider it at a nonperturbative 
level, thus we will assume that there is a result for the flat space amplitude at the nonper- 
turbative level even if we don't know it. Indeed, we are interested in black hole production, 
for which we can calculate the cross section, but not the amplitude, thus in the next section 
we will model the nonperturbative amplitude with a simple black disk that reproduces the 
cross section. Second, the amplitude that we will extract has external legs defined only in 
a flat 4d slice of AdS, not in a general direction in AdS, which is however sufficient for our 
purposes of convolution with ipi according to ()2.6j) . 

Under the assumption that the local string scattering is dominated by the momenta of 
the order of the string scale, l/\foH (which we will shortly see that it's not so innocent as it 
looks), we get by the above ()2.4|) that 



a'p) (2.7) 
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and if \/~&p ^> 1 (high energy scattering in the gauge theory) we see that the integral will 
be concentrated at r scatt ^> r min , where 

V ~ Cf(r/r min )g(n) ~ C(r/r mm )- A g(n) (2.8) 

Then 

A Q cd( P ) ~ I drr 3 (l[( r -^)^)A strmg (pR/r) ~ (-)^'" 4 (2.9) 
J . r p 

as in QCD, and this result was obtained only from conformal invariance. Moreover, for 
states with spin, we replace A with n = Aj — <7, as in QCD. So the scaling with momenta 
comes from the large r asymptotics of the wave function, which itself comes from conformal 
invariance. 

But the last scaling relation treats all momenta the same. As we will mostly be interested 
in the small angle regime, s>t, let's look what happens for the amplitudes as a function 
of s and t. 

Given that A s t r ing — A s t r ing( a 's,a't) we take v = a'\t\ = —a't as integration variable 
(r = u~ 1 ^ 2 r m i n \/&'i) and since s/t = s/t, we have in the new variable (A = >~2i Aj) 

AQCD ~ J dw A/2 - 3 A s tring(vs/\t\,v) 

U°) I Vmligm (2.io) 



_ R W A 4 



Note that the main assumption in deriving this was the large r AdS behaviour of the 
wavefunction (j2.8| h We will have to remember this as a caveat of the formula, but we will 
continue nevertheless to use it, it being the simplest model we can have of an AdS-QCD 
relation. 

In this Polchinski-Strassler form for the QCD amplitude, v max = a't and the amplitude 
to be integrated over is expressed as A(a's, a't). 

Let's observe how a few possible behaviours of the string amplitude translate into QCD 
amplitudes, for future use. We will see in the next subsection that these types of amplitudes 
(power laws plus logarithms) appear from various possible classical cross sections for scat- 
tering in the gravity dual, characterized by b max (s). The various b max (s) of relevance are 
analyzed together in section 5. 

A ~ A 2 s a t l3 log(a't) => 

Aqcd * KA 2 s a t%^) a+ P — * ■ —(log(a't) - 



A/2-2 + a + /T * y ' A/2-2 + a + /T 

fP rJ 1 

A = A 3 lo g (a' S y =, A QCD = KM N2 _ 2 + ^ )Mo 9 (S!s) - A/2 _ 2 + / -ll) 

Giddings [S] points out that as one increases the energy of the gauge theory scattering, 
by ()2.4j) one increases also the relevant energy in string theory. In (j2.4j) . we have seen that 
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corresponding to the string scale 1 / \fa' there is a gauge theory scale 1 / \f& . But there are 
three further (higher) energy scales (in the case when the string coupling g s is small but g s N 
is large). 

The first is the Planck scale 



M P ~ ^ 1/4 (l/v^) = — = N^J— (2.12) 

V 'rain 

Note that we can rescale 4d coordinates such that r m j„ = R, which, since a' = k~ 2 (g s N)~ l l 2 , 
translates also into a! — a'. Then M P = N^A. 

In any case, the Planck scale corresponds in gauge theory to 

M P = g- s llA l\f& = iV 1/4 A (2.13) 

which is the scale at which (real) black holes start to form. 

Giddings [H] proposes that afterwards the black hole production cross section is approx- 
imated by a ~ Tir 2 H ~ E 2 '^ d ~ 3 ' [E 2 !" 1 if we have approximately lOd flat space). As we can 
see, this is based on the simple geometrical picture of a static black hole at a given point in 
AdS, with the cross section equaling its horizon area. We will try to see whether this picture 
is valid. 

The second scale is the string correspondence principle cross-over scale, at which one has 
to stop talking about string intermediate states and instead use black hole virtual states. 
That scale is 



E c ~ gj*/V* = 97 7/4 M P = N 7 / 4 Mp/(g s N) 7 /* = N 2 J^ - ^ 7/d (2.14) 



A 

r~ n {g s N) 7 ' A 



or in the gauge theory 

N 2 A 

E c rsj tL- (2.15) 

The third energy regime is the most interesting, attained when the size of the black hole 
th that was created reaches the AdS size R: 

E ~ M P r 7 H {Mp 2 r d £ 3 ) -> E R = M%R 7 (2.16) 

or in the gauge theory, 

Er = N 2 A (2.17) 

At that energy, the behaviour of th with E changes from (|2.16J) and then so does the cross- 
section. The proposal of jB] is that after this scale we have the onset of the Froissart bound 
in the gauge theory. 

Let us now note that at least after the onset of the Froissart behaviour we have black 
holes being created with Schwarzschild radius greater than the size of AdS, so one might 
ask how come we are still using the Polchinski-Strassler formula (J2.6)) which implies some 
locality for the scattering in the extra dimensions? In fact, in this regime the PS formula 
becomes less and less relevant, since the effective scattering region (where the PS integral is 
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concentrated) is small and close to the IR cut-off (as we will see in the following), and the 
size of the black hole is larger than it. So in this limit, the scattering inside AdS becomes 
more and more classical (the fluctuations due to the PS integral become less than the size 
of the black hole), and the classical cross section calculations in this paper become more 
directly relevant. 

2.1 Black disk calculations 

In the following calculations, we will analyze a classical scattering with black hole creation, 
out of which we get a classical value for a b max (s), and correspondingly a cross section for 
black hole creation, a = ^b 2 max . But in order to use the Polchinski-Strassler formalism and 
relate it to QCD, we must find a quantum amplitude generating the same AdS cross section. 
So we will study first how we get a quantum amplitude out of the classical picture. 

Ideally, one should do a quantum calculation for the amplitude of the AdS scattering, 
not just for the cross section (for which we were able to use the formalism of AS shockwave 
scattering described in the next section) but that would require knowledge of the full non- 
perturbative quantum gravity. As we don't have that, we will try to make a model for the 
amplitude that reproduces the classical cross section calculation, a = nb^^s). 

The simplest thing one can do is to create an eikonal that corresponds to a black disk, 
that is 

Re(S(b, s)) = Im(6(b, s)) = 0, b > b max (s); Im(S(b, s)) = oo, b < b max (2.18) 

Then the eikonal amplitude at a fixed flat 4d slice inside AdS that reproduces the classical 
bmax(s) for scattering constrained to lie in the slice in s,t variables is 

i r rbmax(s) p2tt 

-A(s,t) = -i J d 2 be i<!b (e iS - 1) = i J bdb J d6e iqbcosd 

= 2 m b -^-J 1 (V~tb max (s)) (2.19) 

We will take the amplitude ()2.19|) to be the amplitude for scattering inside AdS space, 
parallel to a fixed 4d slice. Here b max (s) will be determined from black hole production in 
the scattering of two Aichelburg-Sexl Shockwaves inside the curved AdS background. 

In general, the notion of scattering inside AdS space and the associated notions of S 
matrix and in and out states are hard (if not impossible) to define. In particular, the 
original definition of the AdS-CFT correspondence states that a set of good quantities that 
can be related to the boundary CFT are the corellators with legs on the AdS boundary, not 
S matrices. There were several attempts to define S matrices in AdS-CFT by taking the 
flat space limit of AdS, e.g. (THIEVED], but none was entirely satisfactory. What we are 
proposing here is less radical. First of all, unlike these other cases, we are not dealing with 
global AdS, and we are not relating in and out states defined at the boundary of AdS (which 
would be in the Poincare patch at r = and r = oo). Rather, we are looking at scattering 
amplitudes that happen in a 4d slice, at fixed r. Second, when we talk about cross sections 
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we refer to 4d quantities, inside the 4d slice. We are not attempting to define a 5d cross 
section for scattering inside AdS. Here the 5th coordinate of AdS, r, has just the usual role 
of energy scale (according to the usual UV-IR relation in AdS-CFT), as can be easily seen 
in (J2IED. 

We note that if y/tb max (s) <C 1, the result becomes ^b 2 max , so if we take the imaginary 
part of the forward (t=0) scattering amplitude we get 

-ImA e iastic(h, h -> fci, h) = &tot(h, h -> anything) = irb 2 max (2.20) 
s 

But we still need to integrate the amplitude over the AdS slice, using the P-S formula. 
We use the Polchinski-Strassler formula (j2.10|) . i.e. 

Upon inserting the eikonal amplitude (|2.19jl into (|2.10|) . we get 

s y/a!2mK r^=&'t ^ yg /— 

Aqcd = \T\Wy^ I dvv 2 6 - ( ^ )Jl( V ( } 

Let us now look particular cases that will be of interest later on. In all the three cases we will 
study, the result for Aqcd will be of the type mass dimension -2 constant (Kai, Ka\ 4 , 
Kaf) x function of (a's,a't, A and c). Here c stands for a dimensionless number calculable 
in the gravity dual (ai/a', a^ l+2 ^ j 'a' ', a\/a'). Here 01,02,03 are calculable in the gravity 
dual and K encodes the details of the transition to QCD, i.e it depends heavily on the details 
of the IR cut-off of the gravity dual, as can be seen from its expression in (J2.1(Jj) . 
Approximately 4d case: b max (s) = ai^/s. Using 

7M) = r,V lMtfa= a ^(l + V2,2,2 + V2;- a 74) 



4 + 2b 



we get 



2ma 2 K a' 2 s 2 „ .A „ A + 4 a 2 a' 2 st, „ „ ,a 1 * ~ . , 

Aqcd = ~r~ rifj ji 2, — — , —^-) = Ka x x fct. ^, a's, a% A 2.23 

A a' 4 4 4a' a' 

To evaluate the hypergeometric function at large argument (variable), we use that 



Ji(ax) ~ -\l cos(ax + 7r/4) (2.24) 

7TOX 



and then 



J(o,o) ~ * , 1 - [iFi(6 + 1/2; & + 3/2; m)(l + i) - x F x (b + 1/2; b + 3/2; -io)(l - i)\ 
\/ixa lb + 1 



7TO 26+1 
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We get 



, ~/ ,&'t, A-2 , S, 8-A ,Cj' , A-2 8-A _ A „ . 

^ ikat{— ) 2 (-) * = ik{— ) 2 s i t ±a 
a t a 

k = 2v / 2^r(— — )a x 2 cos(^- '-) (2.26) 



2 2 

in the limit ^a' 2 st ^> 1. In the opposite limit, ^a' 2 st <C 1, we get 

2i\ia\K a' 2 s 2 , n 
•Aqcd * — (2.27) 

i 

F/at d dimensional case: b max {s) = a 2 s 2(d " 3) = a 2 s^ 
Using the same formulas, we get 

27rza|Rr a's w A -4 A-4 ^ a't a's 2/3 

^ ra -4/3 + A-2 ( ^ ) aSlF2(1 + 2WTT)' 2 ' 2+ 2WTT)'~T^ ( ^ ) } 

2 

= Ka 2 1+2/5 x , a's, a't, A) (2.28) 

a' 

It has the limit 

- . Pi's ■ <3(A-4) Q;'/; 2-A-4/3 

A?cz? ^ A (— r) 2,3+1 - w) « « 
a a' 

/? = ^ a 2 ^ 2 r( ^ +l) c°s[f(^)]. «i^(^) 2ff »l (2.29) 
(J + 1/2 2p+i 2 2/5 + 1 a' a' 



If d = 5, so (3 — 1/4, we get 

Aqcd ~ fc's"^i^as (2.30) 

Another limit is 

-^i^ri^^ Vs - ^> s<<1 < 2 - 31 > 

Log behaviour: b max (s) = log(s), that would correspond to the onset of the Froissart 
bound, at least in the calculation in [S]. We get 

i — . , .2niaiK f 1 a-s r .„ 
•Aqcd = Ma's) ' / dyy— [Iny + A] J x [B^y(lny + A)} 
\J a't Jo 

a 2 - 

= Ka\ x jctX— n a's, a't, A) 
a 

A = lnfy, B = a 3X m (2.32) 
a' V a 

that unfortunately cannot be algebraically solved exactly. 
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But in the case of A large, c = AB large, we can do the integral, using the large argument 
expansion of the Bessel function, and using 



f 1 , h, / x 1 f 1 , h_i, \ T(b)V2 ,7rb 7i , , 

/ dxx ma; (cos cx — sin cx) ~ — / axa; (smcx— coscx) ~ — cosl ) (2.33) 

and we get 

(a 3 .f^ln(^)y- A _ cos(^(A - 5^ 

* 4v^^r(A - 9/2) 3V ^ [(A - 9/2) cos(^ ( A - 3)) + 2 ^1) ~ 1 

° (2.34) 

Also for A ^> 1, but c = 1, we can expand the Bessel function for small argument 

and do the integration to obtain 

A QCD * a's2malK[^^- - — ^M^)] (2-35) 

Note that for all 3 cases, when we have used the small argument expansion of the Bessel 
function, we have integrated 7rb max (s) 2 with the PS formula. Correspondingly, in the last 
case we studied, we obtain a ~ a|(/n(s)) 2 , as expected. 

We should also note that within this section we have used the simplest model of black 
disk for turning the classical scattering into a quantum scattering amplitude, but for the 
last behaviour (the Froissart bound) there are other choices. For instance, in there are 
given a few forms for the eikonal which when integrated give the Froissart behaviour for the 
cross section. One of them is of possible relevance here. The eikonal 

7rA b^ 
Im 5(6, s) = -^s A exp{-——} (2.36) 
B(s) B{s) 

contains a parameter analogous to the classical maximum impact parameter b max (s) we have 
used: 

B{s) = 2alns + c (2.37) 

where however a and d — A + 1 are related to the Pomeron trajectory ap(t) = d + a t. 
Integrating this eikonal also gives the Froissart behaviour for the total cross section 

a t (s) ~ 4vr a Aln 2 s + o{ln s) (2.38) 

We will restrict ourselves to the simplest black disk eikonal as we did in this section, as there 
is no physical argument on why to choose a particular nontrivial eikonal (in the full nonper- 
turbative quantum gravity one should be able to calculate the actual nontrivial eikonal, of 
course, but in practice there seems to be no reason to prefer a nontrivial phenomenological 
eikonal like the one in (J2.36)) over another). 
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3 Black hole creation in high energy scattering; gen- 
eral formalism 



As we saw in the previous section, at sufficiently high energies in the gauge theory, in the 
gravity dual we will have an inelastic scattering with black hole creation. So we have to be 
able to analyze the black hole creation in a general background. 

In [03] we have extended the formalism in ^l] for calculating black hole creation cross 
sections via analyzing the scattering of two Aichelburg-Sexl Shockwaves inside a curved 
background of AdS type. Here we review the formalism in order to apply it in the next 
section. 

The Aichelburg-Sexl shockwave [221 

d-2 

ds 2 = -dx + dx~ + (dx+) 2 5(x+)$(x*) + ^(cfei) 2 (3.1) 

1=1 

has as a source a massless particle of momentum p ( "photon" ) , with 

T ++ = ^V)<5(x + ) (3.2) 
In flat space, the Einstein equation R ++ = 8nGT ++ implies 

d 2 <S>(x l ) = -l&irGpS^ix*) (3.3) 

($ is harmonic with source). 

't Hooft [23 121] has argued that one can describe the scattering of two massless particles 
at energies close to (but under) the Planck scale, (m^ <C Mp,Gs ~ 1, yet Gs < 1) 

(2) 

as follows. Particle two creates a massless shockwave of momentum and particle one 
follows a massless geodesic in that metric. He has shown that the S matrix corresponds to 
a gravitational Rutherford scattering (single graviton being exchanged). 

At higher energies (Gs 3> 1), one has to consider that both massless particles create A-S 
Shockwaves, and this nonlinear process is hard to compute. At most one can compute the 
metric perturbation away from the interaction point as in j2H]- But one can use a formalism 
to give a lower bound on the size of the black hole being created, and estimate the maximum 
impact parameter that forms a black hole. 

We will be interested in scattering that occurs in a gravitational background of AdS type, 
maybe with a string-corrected source. So we will use a general $ and general dimensionality, 
and a background of the type (the notation is for the one brane RS model, but it is easy to 
generalize for the AdS and 2-brane RS model cases) 

ds 2 = e ~ 2lyl/l [-dudv + dx 2 } + dy 2 (3.4) 

Let us denote e _2 ' y '^ = A and let = Ag^ for the metric in both x and y coordinates 
(transverse). 

We will analyze the collision of two A-S waves in this background, one moving in the 
x + (u), one moving in the x~ (v) direction. The metric in the collision region cannot be 
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calculated exactly even in flat space, but there is an alternative for checking for the presence 
of a black hole in the future of the collision. 

Due to a suggestion made originally by Penrose, one can find a trapped surface at the 
interaction point u=v=0, that is, a closed D-2 dimensional surface the outer normals of 
which (in both future-oriented directions) have zero convergence. By a GR theorem, we 
know that there will be a horizon forming outside the trapped surface, therefore of area at 
least as big as the trapped surface area. 

The metric with a single A-S wave moving in the v direction (at u = 0) in the given 
background is (see [THj for more details) 

ds 2 = e~ m/ \-dudv + dx 2 + $(f , y)5(u)du 2 ) + df (3.5) 

It is useful to perform a transformation of coordinates to eliminate the delta function singu- 
larity in the metric. One finds the transformation 



u = u 

v = v + $0(u) + ^p(<9;$6^ j + A{d y t>) 2 ) 

uQiu 
~2 



x l = x i H ^-g ij dj^ 

y = y + ^Ady® (3.6) 



giving 



ds 2 = A[— dudv + dx 2 + uQ(u)didj$dx l 'dxP\ 

+dy 2 [l + u6(u)Ad 2 ^} + dydx'vBMAdidyQ + dydAu9{u)d y <$> + o{u 2 ) (3.7) 



where 



A = e- 2 ^' 1 + o(u) 2 => A\ U=Q = e~ 2 ^ 1 ; dA\ u=0 = -jA[dy + ^d y $du] (3.8) 

When taking two A-S waves, one in u and one in v, the metric for u < 0, v < (before 
the collision) contains the linear superposition of the 2 waves, i.e. §\5{u)du 2 + ^2^(v)dv 2 , 
and a trapped surface will form at u = v = 0. The trapped surface in the new coordinates 
(u, v) is composed of two "disks" Si and S2 glued on to their common boundary C, namely 

"disk" 1- {v = u = 0}, (*! = on C) 

"disk" 2- {u = -* 2 (f), v = 0} (* 2 = on C). 

The null geodesies through the first disk, {v = — ^f(x), u = 0} are defined by the tangent 
vector 

^ = u-^ + v^+±% (3.9) 
One finds that (by calculating ii, v, x and then lowering indices) 



f = --giidiVdjVdu -dv- d i ^dx i (3.10) 
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and therefore the convergence of the null normals defined by £, 9 = g^D^j is (see ^U] for 
more details) 

= -V 2 (^i - $1) (3.11) 
and similarly for the second surface, and 

V 2 = ^S/l + dl-- l sgn{y)d y (3.12) 

Here $i, $2 are the profiles of the two waves, whose centers are separated by the impact 
parameter b = X\ — x 2 - 

At b=0 and in flat D dimensions (/ — > 00) one can choose = $j, and then as we can 
easily see, the trapped surface S = SiUS 2 corresponds as advocated to the interaction point 
u = v = 0. 

The common boundary C is defined then by ^> = $ = c (constant) or rather, one redefines 
\& — *■ $ — c to have \I/ = on the boundary. We will keep the constant for convenience, but 
remember that we should subtract it at the end. The continuity condition for the normal 
geodesies £ along C gives then 

g ij d^d^ = 4 (W) 2 + A(d y y) 2 = 4 (3.13) 

At b=0 and in flat D dimensions, the \& = $ = c and the continuity condition are 
compatible with the boundary C being a circle (r=const), and then the continuity condition 
fixes the radius of the circle. 

In curved background we can't choose ^ = $ anymore, and the shape of C is fixed by 
requiring that the two equations are compatible. 

Therefore we write for the trapped surface (surface of zero convergence, 8 = 0) 

^ = $ + (^V 2 C = (3.14) 
Now the trapped surface f(r, y) = is defined by both \I/ = C (const) and by g^d^d^ = 

4. 

We will find the shape of the surface perturbatively in y, away from the flat 4 dimensions 
at y=0. Expand ( near y=0 as 

C = Co(r) + Ci(r)y+yC 2 (r) (3.15) 

Then at y=0 V 2 C = implies 

^Co(r)+C 2 (r)-yCi(r) = (3.16) 
and we don't want to upset the flat space solution, so we will take (q = 0. Then 

y = f + ay+£g+... (3.17) 
where / = $^=0, a = d y $\ y=0 + Ci and g = d%$\ y=0 + f Ci- 
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And one has to match the equations 



y = C = f + ay+^g+... (3.18) 

and 

4 = (/' + ya' + |V + ...) 2 + (1 - f + 2^ + ...)(a + yg + ...) 2 (3.19) 

At nonzero impact parameter b, there are two distinct $! and $2 (with different centers), 
and a single curve C, so we can't take = <3>; even in flat D dimensions. The problem 
of finding C together with the functions is complicated, but we have found instead an 
approximation scheme. The continuity condition is now V\l/i • V\&2 = 4, but we approximate 
the size of C by putting ^ = ty(x — Xj), with \1> being the b=0 value. This gives a continuity 
equation 

4 = (1 - ^)[{f + ya' + |V + ...) 2 + (1 - ^ + 2^ + ...)(a + ^ + ...) 2 ]| r=Pc (3.20) 
Now one can say that the area of the real trapped surface satisfies 



S>°\jP 2 c-J (3.21) 
and more importantly we can estimate a maximum b for which a trapped surface forms. 



4 Black hole creation inside AdS and on the RS IR 
brane 

In this section we will apply the formalism of the previous section to the case of A-S shock- 
waves inside AdS and on the IR brane in the 2-brane RS model. The details of the calculation 
of the trapped surface are found in the Appendix, so here we will only show the general fea- 
tures. 

Note that the application of the general formalism for scattering on the UV brane of 
the RS model (as well as for large extra dimensions), and thus for possible applications to 
accelerators was done in [TOj. Here we concentrate on the cases relevant for QCD, namely 
(as we will discuss in detail in the next section) scattering inside AdS and on the IR brane. 

We have already derived the form of the A-S wave inside AdS in [TU]. The function $ 
for an A-S wave centered around r = 0, y = yo is 



Hr,y) = ^±^p e2ie Jl V ° rdqq^J^(qr)K i/2 (^q)I i/2 {^lq) y > y 
(2ir) 2 r 2 Jo 2 

Tdj Ap e2l6 Z / dqq^J^(qr)I d/2 (e^lq)K d/2 (ey^l q ) y < y 
12 r 2 Jo 2 

(4.i; 



(2tt) j 
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Let us therefore first derive the A-S wave solution inside the 2-brane RS metric (the one 
used by Giddings in his calculation of the Froissart bound), with the A-S situated on the IR 
brane. It will be different from the Emparan [26J solution for the A-S wave on the brane in 
the 1 brane RS. We will use the same formalism used for the 1-brane RS case by Emparan 
and for AdS by [HJ]. We would still use the AdS metric 

ds 2 = e~ 2y/l dx 2 + dy 2 (4.2) 

as in [B], but then y £ (— oo,0) and the IR brane is located at y=0 and the UV brane at 
y = —\yuv\- The metric would be valid only in (— \yuv I? 0). For a general metric satisfying 
the required boundary conditions we will take 

ds 2 = e 2 ^ l dx 2 + dy 2 (4.3) 

It matches with the above for its domain of validity (y negative). 

In complete analogy with the calculation in j2B] and ^Qj, we obtain an equation for 
h(q,y) = J d d ~ 2 xe zq ' x $(x,y) which is 

%, V)" + jsgn(y)h(q, y)' - q 2 e - 2 ^ l h(q, y) = (4.4) 

with solution 

Ae-^(I d/2 oiK d/2 )(e~m q ) (4.5) 

and imposing normalizable behaviour at y — ±oo (at the UV brane, if that is moved to 
infinity) we restrict to Id/2- Then imposing the jump condition at y=0 (IR brane), that is 
putting the source of the wave on the IR brane as in we finally get 



2tt J (2vr)^ 2 qI d /2-x{lq) 



(2vr)^ Jo ~ h/2-i{lq) 

(in d=4, AGd+iP = R s l)- By comparison, Giddings has an h 0Q (Newton potential) obtained 
also from a wave equation with sources on the IR brane [B], except his source was a static 
mass (black hole), whereas for us it is a photon 

and since J v {ix) = I v {x) and y is negative, we have almost the same solution, except for us 
the integration is only over d-2 transverse coordinates (as for a massless particle), whereas 
for Giddings the integration is over d-1 transverse coordinates, as appropriate for a massive 
particle. 
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Let us then analyze the AdS scattering. We are interested in the large r behaviour 
[e y l l ljr <C 1 and y = y^) which is relevant for the Froissart bound that we are after. In that 
case, 

$ = ^! e f (2y+yo) ( 48 ) 
We find (see the Appendix for details; here e = y — yo) 

f = 4>|„„ = C^e 6 *"" 

9 = 3?4|,.o + jCi = + J Cl 

As we have explained in the Appendix, we actually get two solutions for the trapped surface. 
There is one solution for which a is neglijible, thus that solution corresponds to what we 
would have if the scattering was four dimensional (and just $ was obtained from the 5d 
equations). But in the r 3> / limit, the AdS warping is very large, and energy scales become 
larger away from the 4d slice, thus we expect the size of the black hole being created (and 
thus of the trapped surface) is increased. And so we will take the solution that takes into 
account the 5d scattering, solution that implies a larger horizon. 

The same situation will be encountered in the case of scattering on the IR brane. The 
space will then be very non-4-dimensional in the r 3> I limit (in a sense, it will be anti- 
4-dimensional), thus we expect that the size of the trapped surface will be increased also 
with respect to pure 4d scattering. This situation is to be contrasted with the scattering 
on the UV brane we analyzed in [HI], m which case for r ^> / the space was approximately 
4 dimensional (the warping is going down away from the brane), and correspondingly we 
found a solution which had just small corrections to the four dimensional scattering. 

The continuity condition for the larger trapped surface is thus 

\ a \f e -yo/i = j ( 4 10 ) 
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so that 



le^&]^ (4.11) 



l l e yo/i 

and since in 4 dimensions R s = 2Giy/s, we have that r max ~ s 1 / 12 . 

If we introduce a nonzero impact parameter the continuity condition becomes 

and with r 2 = x, (3R s l 5 e 5yo ^ 1 ) 2 = a we get the equation 

b 2 

x 7 -ax + a— = (4.13) 
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The maximum b for which it has a solution is (using f(xo) = for f'(xo) = 0) 

2 = 12ra 1/6 1/12 3Rs_ 1/6 _J_ /12 f414 x 

so that 6 maa; ~ as 1 / 12 as well. 

We will now analyze the scattering occuring on the IR brane, the details of which are in 
the Appendix. The wave profile is 

$(r,y) = RJe-W f°° dqJ (qr) Hq) (4.15) 

Jo h{lq) 

and we find that we can use the contour integration over the complex q plane to calculate 
/ = $(r; y = 0) * ^Tjg^ ^^(ii.n) ^ ^ R JMc^ (4.16) 

where g n = ji, n /l are poles of the integral in the complex momentum plane given by the 
zeroes of the Bessel function J\, ji >n . In (I3.18j) from the general formalism, i.e. 

* = f + ay + g^ + ... = C (4.17) 

we have (notice the sign difference in g, due to the IR brane metric ()4.3|) vs. the UV brane 
metric ()3.4j) used in section 3) 

4 4a 

a = Ci, <7 = ^$|y=o-yCi,= -/9?- T (4-18) 

Again, as in the AdS case, we find two solutions, one that would be there if we had a 4d 
scattering, and a larger trapped surface that aapears only when we have 5d scattering. In the 
Appendix we treat in detail the case of the 4d scattering, for completeness. Taking instead 
the larger trapped surface, 

Jr 2 



a = -3^— (4.19) 
r 



we get the continuity condition at y=0 



v 2 ' q x l 2 2 V r q\l I 

It has a solution that is approximately 

1 / ( A \ 
ru — ln[ 



A=( ^ 1)2 T ( i L)2 (4 - 21) 
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We can see that there is a solution by considering the function 

~ g {r) = r - Ar 2 e-^ ir => g'(r) = 1 + (ft - -)2Ar 2 e~ 2 ~ qir (4.22) 

r 

and the solution for r# is given by g(r) = 0. Since g'(r) > if r > 1/ft and 

, 1, 1. R 2 18nC 2 , 
~g(-) = -[l--^-—-^}<0 4.23 

for sufficiently large R s , there will be a solution. 

At nonzero impact parameter, we get the equation 

Its solution is the zero of the function 

ft 2 

^(r) = r - (1 - — )Ar 2 e~ 2qir (4.25) 

And now we have an analysis that is a bit more involved. 

b 2 

g'{r) = 1 + 2ft (r 2 - —)Ae~ 2 ^ r - 2rAe~ 2q ~ ir (4.26) 

We are however only interested in the maximum value b max for which g(r) = has a solution. 

If r 2 < b 2 /2, then #(r) > 0, and g(0) = b 2 A/2, g(b/y/2) = b/y/2. To have a solution of 
gir) = 0, we need a minimum, g'{ri) = 0, with g'(ri) < 0, so necessarily r\ > b 2 /2. 

But if 6 2 /2 <C r 2 , the second term in g'(r) is larger than the third (as ft ~ 1/7 3> r), so 
g'(r) > 0. Thus we need instead 

_6 2 
2r 2 

so if p(ri) < (but close to zero, so that the 1 is negligible in g'(r) and we can use g'( r i) = 0); 
we get 

g( ri ) ~r!(l--e- 25ir ) (4.28) 
ft 

Atb = b max , g(n) = 0, so 

ri ~ — Zn(— ),6 2 ~ 2r 2 
2ft ft 

£w(s) ~ v^Zn(^) = ^/ n [fl,ft(4^)] (4.29) 
2ft ft ft v^Jii 

where i? s = y/sG^, G^ = l/(ZMj> 5 ). Thus we get the same formula for r^(m) as [0] (modulo 
different constants), and thus the classical cross section in the gravity dual for scattering on 
the IR brane is of the expected ln 2 s functional form 

o- = tt&LJs) ^ 27r[l/n(K lv ^ftG 4 )] 2 (4.30) 
ft 



- 1 - e ~ 1 (4.27) 



18 



5 QCD scattering regimes; the Froissart bound 



Finally, now that we have done all the calculations needed for the scattering in the AdS 
dual, let's put everything together. In the gravity dual, the scattering happens inside AdS 
or, in the Froissart regime, effectively on the IR brane. 

In JTU], we have analyzed 't Hooft scattering in AdS, and we have found that we can 
calculate analytically the scattering amplitude in two limits: 

A AdS - ^ * ^fV^/^o/ i); y ^ < j, r « j, g 4 s « 1 

A AdS oc G 4 l 6 st 2 ln(t), r > I (5.1) 

to be compared with the result in flat D dimensions (which can be used in the case y = 
yo, r <C I, for instance) 

A AdS ~ -j- (5.2) 
All cases still give A ~ G^s, only the t behaviour is modified, so 

da Ads ^da A ds ln -.2 / e o\ 

= S^n~ ~ (G4S) (5 - 3) 

Since we have a well defined amplitude, we can use the Polchinski-Strassler formula (|2.10j) 
to relate to a QCD amplitude. Of course we have to remember that the formula is valid only 
for the large r region, away from the IR cut-off, where we can approximate the metric with 
AdS and the wavefunction dependence of r with the power law r _A . But unfortunately, in 
all the cases of interest, studied in (J2.11j) . the main contribution to the integral comes from 
the region of large v {y ~ u m ax = &'t), corresponding to the region of small r, r ~ r min . 

This is however not so bad as it seems, since this IR modification will translate only in 
the modification of the factorized t behaviour (coming from the integration over u), while 
the dependence of s comes from the s dependence of the AdS amplitude. We can check this 
by looking at the formula (|2.1()j) . As we will be interested only in the leading s behaviour 
and not the t dependence (be it multiplicative or not), our results will still be valid (see 
also the discussion further of the original Polchinski-Strassler paper [3] later on, especially 
around (|5.11jl ). 

We also see that in all the cases of interest treated in (|2.11j) . the effect of the AdS 
integration is only to change the overall normalization, as well as the subleading behaviour. 
The leading behaviour is kept the same. 

In conclusion, for the case at hand, of the 't Hooft scattering in AdS, where A A ds ~ G4S, 
by using the PS formula (J2.10J) we will still have the same s dependence 

A QCD ~G A s^ d -^ = \^ ~ \A QCD ? - (G 4 sf (5.4) 

As we mentioned, the relation between AdS and QCD energy scales is: y/a's\ A ds < 
V&s\qcd and \fc/t\ AdS < *f&i\QCD- Here 

« = A2 (5-5) 
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where A is the mass gap (the lightest glueball state). And the minimum r in the cut-off AdS 
is r min = R 2 A, R = I in our notation. 

We can rescale the 4d coordinates x such as r m j„ = R =>- a' — a'. The Planck scales are 

-1/4 -1/4 

M P (AdS) = ^M P = = N 1/4 A 
vat' yd' 




J-<£- G^Mf = -L- (5.6) 
Mp~ M P P VNA 2 

So from 't Hooft scattering in AdS we get a Rutherford-type behaviour in QCD as well, 
except with the effective coupling G 4 s. 

This behaviour is universal, and if we interpret it as single particle exchange between the 
gauge invariant Yang-Mills states (glueballs), we see that the particle being exchanged must 
also be colourless. 

So in the energy regime G^s < 1, but at energies larger than the mass gap, i.e. for 

= AN 1 / 4 (5.7) 



colourless states should obey Rutherford scattering behaviour, with a universal colourless 
single-particle exchange, with universal effective coupling G^s. In the case of real QCD 
(N=3), the energy regime is nonexistent, thus we can't draw any lessons from experiments, 
and the interest in this regime is therefore just theoretical. 

However, as the universal coupling was obtained from a spin 2 exchange in AdS (graviton), 
it is natural to assume that the same thing happens in the gauge theory, namely a universal 
spin 2 colourless particle is exchanged in this regime that behaves as a graviton. There 
is a natural candidate for such a composite "particle", namely the gauge theory dual to 
the graviton, the energy-momentum tensor. It is not obvious why in this energy regime 
the energy-momentum tensor should have a graviton-like universal coupling to all colourless 
states, nor why graviton exchange in the bulk should be dual to energy-momentum tensor 
exchange on the boundary, but this is the only plausible candidate. 

This is an elastic scattering in AdS, and it is an elastic cross section in QCD. It could 
then maybe be possible to detect this elastic a even at higher energies, when the amplitude 
is mostly inelastic. 

As we go even higher in energies {pis ^> 1), we will start observing the Regge behaviour 
noted by Polchinski and Strassler Then in AdS we need to use the string Virasoro-Shapiro 
amplitude for massless external states (s+t+u=^ mf = 0) 

n r(-q'x/4) r- r(aW4 + aV4) 

iL r ( 1 + «W r(i + ay4) 

x - T ^ 9 ^ , - r y^K( P ^) (5.8) 



T(l - a's/4 - a't/4) T(l + a't/A) 
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becoming in the small angle (s>t) region of interest 



^aV^-H-^r'^jQ^Apol. tense.) ~ (a^I^lB. 

(5.9) 

Regge behaviour in QCD was found to correspond to Regge behaviour in AdS (to the flat 
space Regge behaviour of this VS amplitude, to be exact). The approximation used though 
is the same approximation that we needed to use for the 't Hooft scattering. Namely, in 
general the v integral is dominated by (using stationary phase approximation) 

u ° = ~WiW) (5 ' 10) 

But Regge behaviour is obtained in the case that is > u max = a't, when the integration is 
dominated by the upper limit of the integral, corresponding to r close to r min . 

But as in the case of 't Hooft scattering, we only need to assume that the maximum of 
the integral is still outside the region of integration, so that we can approximate the integral 
with its upper limit. Presumably the existence of the mass gap is enough to satisfy this 
requirement. Then the factorized t dependence will be modified, but the s behaviour will 
still be of Regge type: 

A ~ (a's) 2+& ' t/2 (5.11) 

Finally, for the elastic amplitude, [3] find yet another regime. If v max > u (to be rigorous 
we would need v max ^> z/ , corresponding to r scaU 3> r min ), that is if 

(A -4) , 

a't > ; . T . . 5.12 

ln(s/\t\) K ' 

one obtains 

^~s 2 |t|- A / 2 [/n(5/|*|)] 1-A/2 (5-13) 

Let us now turn to energies higher than the Planck scale in AdS, s > Mp s > Mp, 
when we will produce black holes in the scattering. There are 3 dimensionless parameters 
characterizing the scattering in AdS, R s /l, G^s and /Mp 5 , so let's first derive their relations 
to QCD variables. 

Since Mp j5 = N l / A /I in our notation, then 



so that 



I y/N M P 

Also, since G 4 S = R\jG^ we get 



^<4v(#)lQ^ (5-15) 



G A S < iV" 1/4 (#) 2 |QCD (5.16) 



'Mr 



21 



Equivalently, in terms of &', the two parameters satisfy 



R ^ 4 ^ 2 

T ~ 7^ v/ ^' QCD ' G * S ~ J^TJl( & ' s ^QCD (5-17) 



Finally, /Mp j5 = iV 1 / 4 ^> 1, so one first reaches the AdS scale l/l, then Mp = Mp j5 , and 
then the scale Er = Mp(lMp) d ~ 3 , when the black hole size is comparable with the AdS size. 

Note that the dimensionality of the gravity dual plays an important role. Most of the 
calculations in this paper were done assuming that there is only AdS?, and forgetting about 
the compact space. In the case of 't Hooft scattering, we argued that the relevant behaviour 
was independent of the dimensionality of the space. In the case of black hole creation, the 
dependence is more important. 

Before the size of AdS becomes important, black hole formation can be approximated as 
being in flat space. In that case, we have calculated in a lower limit on the maximum 
impact parameter that creates a black hole in d dimensions, 

, < 2[ - ,^ = r W {D _ 3) 



max--i D _ 2 l D 



[D - 2]^=3 



' L (D-2)fW 



4fi 



D-3 



(5.1J 



In any case, we see that b m£kX (s) ~ as 2 '- - 3 ') = as 13 (a= constant). In section 2.1 we have 
used a simple black disk model to create an imaginary elastic scattering amplitude that was 
substituted in the Polchinski-Strassler formula (j2.1()j) . We have derived the forward (t=0) 
imaginary part of the amplitude, giving us the total QCD scattering cross-section in this 
regime 

na 2 K .a's.na na 2 K . 61 s . i ,„ 

@qcd = 2/3 + A/2 _ 1 (^7) = i /(D -3) + A/2-l ( ^ )g " 3 (5 ' 19) 

As we have explained in section 2.1, this is of the type 

2 

2 - 

(KawP) x fct.( , a's, a% A) (5.20) 

a' 

and here a 2 ^ 1+2 ^ /a' is a number depending on the gravity dual (dimension D, Newton 
constant G), whereas Ka 2 ^ 1+2 ^ is a dimension -2 constant that also depends strongly on 
the details of the IR cut-off (through K). 

We see that the higher the dimensionality, the smaller the dependence on s. For d=5 we 
have a ~ s 1 / 2 , whereas for d=10 we have a ~ s 1 ^ 7 . 

At even higher energies (y/s > Er), we will start feeling the effects of the AdS size. We 
have calculated in section 4 that the maximum impact parameter for black hole formation 
is at least equal to 

h = T 1,12 xr^lc yoll \ Rs 1 1/6 = - \l^r (5 21) 

Umax ' V 7 [ ley°/ 1 71/1222/3 Y 7 max 1 } 
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where R s = 2Giy/s, so that b max — a's 1 ^ 12 . 

Then using the same black disk model to substitute in the P-S formula, we get 

txo! 2 K a's 2(3 txo! 2 K a's.i . . 

aQCD = 2/5 + A/2-l ( ^ ) = l/6 + A/2-l ( ^ )6 (5 ' 22) 

But as in the previous cases, the P-S formula shows that the main integration region is 
near the cut-off r min , which corresponds to the IR brane in the 2-brane RS model. But in 
that case we can't approximate the scattering as being in AdS, since the 4d size of the black 
hole formed is comparable to the AdS scale. As p ^> / the horizon will stretch over a size 
Ay > I, and if y ~ ym it would look as if the black hole is approximately on the IR brane. 
So the approximation of A-S in AdS will break down and instead the good approximation 
would be the two A-S Shockwaves being on the IR brane. 

So it is not even clear that there is an intermediate regime of the type in (|5.22|) . but it is 
clear that the cross section will begin flattening out, finally to settle into the final behaviour, 
corresponding to scattering on the IR brane. 

In the second part in section 4, we have calculated that for A-S scattering on the IR 
brane, we get a maximum impact parameter for black hole formation that is at least equal 
to b max = y/2ri, with 

r 1= -U4); -4=(^)^& (5.23) 
2<?i qi 2 / qil 

In that calculation we have used the metric 

ds 2 = e 2 ^ l dx 2 + dy 2 = ^(-^dx 2 ) + -dr 2 (5.24) 

1 ' min ' 

so to go back to the real coordinates we substitute 

I piused) 
pyreal) = p(used) = — — — (5.25) 

Thus 

bmax{s) ^ ^-V2^-ln(— ) = —InlRsq^ 3 /^ /2 )} 
IA 2q x qi q x \/2ji'{ 



(5.26) 



and R s = ^G^, G^ = l/(ZMj> 5 ). As advocated, we get the same formula for r^(m) as in 
[Oj, modulo different constants. 

The gravitational cross section is 

a ~ 7r [^i/n(K v ^g 1 G 4 )] 2 (5.27) 

and ctqcd ~ as we argued. 

Thus as expected, the final behaviour of the QCD scattering amplitude corresponds to 
the Froissart unitarity bound. We have obtained the same behaviour that Giddings [H] has 
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proposed, but in a more rigorous setting. Let us compare to the calculation in 0. There, 
the "Newton potential" h 00 was calculated in linearized gravity , obtaining 



h °°~ e2 'J Wy^ e pj^m (5 ' 28) 

which was then used for an estimate of the horizon size by h 00 ~ 1 =>- r = r#, and a 
geometric cross-section approximation a ~ vrr^ was used for the black hole. 

In our case, we obtain the exact A-S shockwave solution on the IR brane, which is 

*(r,y) = ^f^ f ^ 



2tt 7 (2vr) d - 2 qI d /2-i(k) 



(27TJ 2 r 2 Jo -'(i/2-1 

and since J u (ix) = I u (x) and y is negative, is very similar to the Giddings case, except 
that for us the integration is only over d-2 transverse coordinates (as for a massless particle), 
whereas for jH] the integration is over d-1 transverse coordinates, as appropriate for a massive 
particle. The /loo in [S] was obtained also from a wave equation with sources on the IR brane, 
except there the source was a static mass (black hole), whereas for us it is a photon. 
At y=0, one obtains similar behaviours, 

h 00 ^—Y e-^ J J^lA „ 1 -ftr^fc) (5 . 30) 
zTir z — ' lcl\ n 2-nr la\ i 

versus 

$(r,y = 0) ~ RsJ^Y, • 7l "" /2j2 °' 1 ' w) e-^ ~ R^^C^ (5.31) 

n ' 

and in both cases the fact that allows the logarithmic behaviour of r# is the exponential 
e~ ?ir , itself coming from the presence of the pole in the momentum space integrand. In our 
case, we have the advantage of the scattering picture, in which we can calculate directly the 
cross section for black hole creation. 

Finally, in ^U] we have addressed the issue of string corrections to the scattering. This 
is very relevant for the case of QCD, since string a' and g s corrections in AdS translate into 
1/N and 1/ '{qymN) corrections in the gauge theory, bringing us closer to the case of QCD. 
So it is important to realize their effect. 

In [TO], string corrections were analyzed using a formalism of Amati and Klimcik [27| (as 
well as using the formalism in j2H], based on the action in but we will not describe it), in 
which one obtains a string-corrected A-S metric. These corrections depend the dimensionless 
ratio (in d=4; here Y = a'logipt's)) 

R2 s 9 2 9 2 a's . . 

Y log{a's) (4vr) 2 { 1 



/ dqq 2 J^( q r)— — — 5.29 

Jo 2 I d /2-i{lq) 
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When this parameter is large, the shockwave is approximately A-S, with exponentially 
small corrections, namely 



m = -^V(2 log A_ e -^ ( ^)-i + ...) = - Rs( 2 l 0g ±.-e-&(^ + ...) (5.33) 

and the maximum impact parameter for black hole creation also increases, but with expo- 
nentially small corrections: 

B max = ^=(1 + e~ R2 ^) (5.34) 
When R 2 S /Y <C 1, the shockwave is not of A-S type 

m = -2*.(j^4 - WJW^j + -) (5-35) 

and we can show that B max ~ -v/F ^ R s , (but we can't compute the actual value), so one 
has a huge increase in the cross section. 

So string corrections can be quite important, but in the case of s — > oo (the Froissart 
unitarity bound), R 2 S /Y ^> 1, and the corrections discussed here are exponentially small. 
Thus the Froissart unitarity bound is unaffected, as expected. 

While this calculation was in flat d=4, not in the gravity dual, and string corrections 
were described using a very simple model, the smallness of the corrections, together with the 
fact that we do obtain what we expect, namely the Froissart bound, makes us believe that 
we are correctly describing a QCD phenomenon. 



6 Conclusions 

In this paper we have analyzed high energy QCD scattering in the small angle region s> t. 
We have applied the high energy gravitational scattering calculations in ^U] to the simple 
model of QCD gravity dual used in j3J. Namely, for high energy QCD many of the observed 
features can be deduced from just a AdS dual cut-off in the IR (and maybe in the UV), 
giving a 2-brane RS model. The gravitational scattering was analyzed using a shockwave 
analysis. At energies smaller than the Planck scale, the scattering of two massless particles is 
described by null geodesies propagating in the shockwave background ('t Hooft scattering). 
At energies higher than the Planck scale, we need to take two A-S Shockwaves, and black 
holes are being formed in the future of the collision, for an impact parameter less than a 

't Hooft scattering corresponds in the gauge theory to a very restrictive regime (A < 
y/s < AN 1 / 4 ) that is too restrictive for real QCD (N=3). For N large though, we obtained a 
Rutherford-type scattering with effective coupling G 4 s, implying a universal single-particle 
exchange. We conjectured that this comes from an exchange of a universal, graviton-like 
spin 2 "particle" being exchanged, most likely the energy momentum-tensor (the dual of the 
graviton). As the cross-section is elastic, one could however maybe detect this elastic a even 
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at higher energies in the gauge theory, when the amplitude is mostly inelastic, thus being of 
possible relevance to real QCD. 

At higher energies, Regge behaviour sets in, as described by [3J. Regge behaviour of the 
flat space Virasoro-Shapiro amplitude translates directly into Regge behaviour in QCD. 

At even higher energies, black hole production sets in inside the gravity dual, and it can 
be described as if happening in an approximately flat d-dimensional space, giving the power 
law behaviour in (|5.19j) . namely oqcd ~ s 1 /^ -3 ). 

There is a possible transition region when the size of AdS becomes important, and the 
scattering in the gravity dual created black holes inside AdS, giving again a power law 
behaviour (|5.22jl . namely <Jqcd ~ s 1 ^ 6 . 

Finally, in the last energy regime most of the scattering happens on the IR brane, and 
the size of the surrounding AdS is important. We obtain a shockwave solution that has the 
same exponential behaviour h ~ e~ 5ir as the linearized black hole solution used by Giddings 
The scattering of two such Shockwaves gives a scattering cross section saturating the 
Froissart bound (|5.27|) . 

We have looked at string corrections to the scattering of two modified Shockwaves, and 
we have found that in the simple Amati-Klimcik model used in ^U|, and in flat d=4, we 
obtain exponentially small corrections in the Froissart limit s — > oo. This makes us believe 
that even in the case of real QCD, when N and gy M iV are finite, thus we have large string 
corrections in the gravity dual, our calculation of the Froissart bound still applies. 

Acknowledgements We would like to thank Matt Strassler for a lot of useful discussions 
on his work, and also thank Antal Jevicki and Carlos Nunez. This research was supported 
in part by DOE grant DE-FE0291ER40688-Task A. 
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Appendix A. Trapped surface calculations 

Let's apply the general formalism to calculate the shape of the trapped surface in AdS. 
The continuity equation at b=0 would be, for ^ = $, 

(O^) 2 + e~ 2y/l (d y $) 2 = 4 (A.l) 

which for e y ^l/r ^> 1 gives 

3y-y C 2 r 1 9 1 . . 

6 ' ^4T \r 2 + l 2 (ev/ 1 - e W*) 2 ] 2 + 4ZW[r 2 + Z 2 (e^ - e^A)2] J = ^ ' ^ 

(where C* = 2R s l 2 ), since 

^ e 2i 1 

* ; ; A -3 

2/ v / r 2 + / 2 (e^ -e^A) 2 

, /? / Q r 2 

1 2r 2j 1 4/ 2 e 2 ^ J 1 1 

where as we can see the second term (coming from (d y §) 2 ) is a small correction. 
For l<rwe have 

$ = ^L e H*y+vo) (A.5) 



At y = y we get 



That implies 

^ e f( 2 ^o) ]2[1 



r 6iU 6 a,,.^,,,. 4 r 2 



r 



7 



9/ 2 e 2 ^ J 



However, now the second term (coming again from (c^$) 2 ) is dominant. 

But we still need to modify $: ^ = $ + £, as in the RS case. Using the formulas from 
the RS case described in section 3, but remembering that now we expand in e = y — yo, not 
in y, we get 

V = C = f + ae+ e -g+... (A.7) 

where 

f = $\ v = V0 , a = ( 1 + d y $\ y=yo (A.8) 

Using 

® = K-^ / dqq^J^{qr)K i/2 {e^lq)I A/2 {e^ l lq) (A.9) 
r 2 Jo 2 

and zK' v (z) + vK v {z) = —zK v _i(z), we get in d=4 

«9 e $ = — / rfggJ (gr)(e e/ 'e yo/i /g)A^ 1 (e e/ 'e J/o/ '/g)/ 2 (e s/o/i /g) (A.10) 
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Then using h(z) = —21\{z)/z + Iq(z) and 



dxx 2 J (x)K l (ax)I (ax) = 2a(l + 4a 2 )~ 3/2 

Jo 

we get 

Ke- 2y °/ 1 6e 4y °/ l l 4 /r 4 + 6e 2y °/ l l 2 /r 2 + 1 
£$ ' e=0 ~ /3 ^ (l +4e 2 ^/ i / 2 /r 2 ) 3 / 2 ^ 

At e y °/ l l/r < 1, using if = Ce 2y °/ 1 = 2R s l 2 e 2y »'\ we have 



/3(f) 

£ $L=o ~4Ce 6yo// — =4-, 



e=0 



Then also 



Q POO 

<9 e 2 <t>| £=0 = ^$| e =o + Ke 2y °/ 1 J dqq 3 J (qr)K (e y °/ l lq)I 2 (e y ^lq) 

Now however we can only do the remaining integral if the argument is small, 
could try using again I^z) = —2I 1 (z)/z + Io(z) and 

3l/u// s T , s , 1 + 2a 2 + 4a 4 
x Jo{x)K {ax)l {ax) = — 4- 



(l + 4a 2 ) 3 / 2 

but the Ii integral cannot be done. Instead, expand K (ax)l2{ax) and find 

POO 

/ x 3 J (x)K (ax)I 2 (ax) ~ 8a 2 - 96a 4 + o(a 6 ) 
Jo 

and thus for e yo ^ l l/r <C 1 we have 



/4 

<9 2 $| e=0 ~ 16 • 2R 8 — e 6yo/l 



rp6 



Thus at e y °/ l l/r <C 1 and y — yo 



4 

a = T / + C, 



If a is nonzero, we have to match 

,2 



4 = /' 2 + aV 2yo// + e(2a7' - 2y e - 2W ' + 2age- 2yo/l ) + ... 
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with 

C = f + ae + ... (A.20) 

If we put (i = we then get 

1 = (jRs i! e %o// )2 ^ e -2W* (1 + 9 il e W + 6j + ...) = {Lf^e-^H^ + 9 il^/i + ^ + 

(A.21) 

where now the first term comes from (<9^\I/) 2 , to be matched with 

C 2 = / 2 (l + y + ...) (A.22) 

which doesn't work. The next try is to put a nonzero (j that still keeps a nonzero. If it 
keeps it of the same order as the first term, we will have 

a = -jf (A.23) 

which implies we have to match 

1 = (|)2^_ e -2yo/i (1 + 6| + ...) (A.24) 

with 

C 2 = / 2 (l-^ + ...) (A.25) 

This happens to have a solution, a = —3! But let us see if it is unique. We could also have 
a = to this order (a = 0). Then let's assume that a is proportional to the next term in 
the expansion, namely 

a = pt e vo/i (A.26) 



r 

so that Ci = 4/// + a and so g = Aa/l + .... Then one finds 

1 = ^ 2 ^ 2+ ^ 1 + iwtw ) (A - 27) 

to be matched with 

C 2 = f 2 (l + 2e^-e yo/l + ...) (A.28) 
r 

and as we see the functional dependence is different. For a higher power in the expansion if 

a = pl e w>/i(L e vo/i)» n > i (A.29) 
r r 

the first equation will be 

1 = (^) 2 {1 + ^-[f3 2 (- e yo/l ) 2n + 2(3(7 + n)(-e yo/l ) n+1 }} (A.30) 
r 61 r r 
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whereas the second will be 

C 2 = f 2 (l + 2e^e yo/l ( l -e yo/l ) n + ...) (A.31) 

We see that for n > 1 we have a mismatch, we would need 2(n+7)/6 = 2, which is impossible, 
whereas for n = there is a different functional dependence (and for n < — 1 we can treat it 
separately and convince ourselves that the functional dependence is also different). But for 
n=l we have another solution, the matching condition is (3{f3 + 4) = 0, so 

a = -4^ (A.32) 

Finally, we can't have a=0 exactly, since then we need to match 

A = r 2 + y 2 {fg' + g 2 e-^l l ) + ... (A.33) 

to 

C = f+£g+... (A.34) 

and then we have g=0. 

So we have two solutions, a = 3/// and a = — 4/7/r 2 . Which should we choose? A 
physical argument shows us what happens. In the second solution, a is negligible at y — yo, 
so we have the same continuity condition that we would have if the shockwave scattering 
was 4 dimensional. But now we are in AdS and the warping is very large, so we are at higher 
energy away from y = yo, therefore we expect that the black hole formed is larger than what 
we would have in the four dimensional case. So whereas the second solution describes a 
trapped surface that would be there even if the space would be four dimensional, the second 
trapped surface is larger and is due to the very large warping outside the 4 dimensional slice. 

We might be worried that there is some theorem stating there should only be a trapped 
surface, as the trapped surface problem is similar to the Green problem with Neumann 
boundary conditions in electrostatics. But this is not quite so, since now we must also 
determine the boundary C from the condition that (V^) 2 = 4 matches \& = (arbitrary!) 
constant, together with the solution to the Laplace equation A(\l/ — <&) =0, so it is not quite 
the same Green problem. 

Thus we take the larger of trapped surfaces, with a = —af/l, in which case the continuity 
condition becomes 

Ulf ., 

e~ yo/l = 1 (A.35) 
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so that 

r = w =/<»/' (A.36) 

and since R s = 2Gi^/s, we have that r max ~ s 1 / 12 . 

Now let us do the same for the case of the wave on the RS IR brane. 
We have 

,_, Jd-i(qr) 

dn d - 3 e^ cose = (2tt)— * (A.37) 

(qr) 2 
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which we can apply for 

J d d - 2 qe iq ^f(q) = J q d ~ 3 dq J dn d _ 3 e iqrcosd /(g) (A.38) 

However, in the case of P, for the calculation of h 00 , the integral that one has corresponds 
formally to d=5 in the above, but it is more useful to do the integral in a different way, 
namely to write 

/°° f —1m f°° 

q 2 dq J dtt 2 e iqrcose f(q) = — — J qdq{e iqT - e- iqr )f{q) (A.39) 

and if f is even, i.e. /(g) = f(—q), we can rewrite it as 



-2ni 



+00 



OG 



r 

For jH], on the IR brane, that is at y=0, the function f is 



qdqf(q)e iqr (A.40) 



_ Uiql) _ h{ql) (A4 ) 

f[q) ~ iqj x {iql) ~ qh(qD ( } 



(I v (iz) = i v J v {—z)). More precisely, 



M , = 0) * /^-- J2(i?i) 



(2ir) 3 J r iqj^iql) 
1 :( Z™) (A.42) 



(2irf r 'i-00 lh(ql) 

But there is a theorem: For a complex function f(z) such that lim^oo f(z) = 0(lm(z) > 
0), and a real a > 0, we have 

/+00 
f{x)e iax dx = 2m ^ Rez[F,a], F{z) = f{z)e iaz (A.43) 

00 lm(a)>0 

In Giddings's case /(g) = I 2 (ql) / h(ql) 
We get 

^^^E 6 "* 1 "— (A.44) 
2vrr ^ ai „ 



it 



where we have defined q = iq and the behaviour of the Bessel function near a pole is 

Jx(z) ~ a lsn (z - z n ), z ^ z n (A.45) 
The zeroes of J\ are called ji >n , so g„ = ji, n /l. Then 

Z7rr z — ' <ai n Lixr law 

n ' ' 

31 



Let us come back to our case, of computing $(r, y). We will first try to do the integral ex- 
actly, and see that unfortunately we get nonsensical results, and then make an approximation 
that allows us to do the integral. 

The integral we have is 



/*OG /*Z7T /*oo r-\-7V / z 

= qdq d9e iqrcos0 f(q) = / qdqf(q) / d6(e iqrco8e - e - iqrcoa6 ) 

Jo JO Jo J-w/2 

/+7r/2 /" + 0O 

dO / qdqf(q)e tqrcos0 (A.47) 
-7r/2 «/— OO 



-tt/2 p+oo 
r/ 

and we see that in the 9 integration regime 1 > cos 9 > 0, and 



as before, so that 



and correspondingly 



qh{ql) 



+W/2 d92m V e-ffnrooee^n) (A 4g) 

tt/2 „ «*l,n 



(2tt) 2 ^ ai,„ J_ x/2 

= T^^E^^^O-^r)) (A.50) 
(271-^ ^— ' ai n 

We are interested in the limit q n r 3> 1, and 

L (*) ~ I (z) ~ -r^e 2 (A.51) 

V27T2 

at large z, so that is not very useful. We can instead expand the integral already and get 

/ + 1 d9e- qnrcose = 2 / + ' d9e~ q ' nrcose ~ 2 / dfcT«" r * ~ — (A. 52) 

J-tt/2 Jo Jo Qn r 

so it would seem that there is no exponential behaviour! However, the approximations used 
were contradictory, since we have expanded the previous integral about the point where 
cos# = 0, which is exactly the point where the contour integration theorem doesn't work. 

So we must find an approximation regime when we can do the integral. 

This time we do the angular integral and obtain 

I = 2n dqUqr) 1 ^- (A.53) 
Jo h{ql) 
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Since we want to have r > we can use the large argument expansion of Jo, 

[Y 1 

Jq(z) ~ \ — cos(z — 7r/4) = (cos z + sinz) (A. 54) 

V 7i z ynz 

thus 

2~ r° , _-, /2 e igr ~ i7T ^ + u q n 

,-rl dqq -' 2 W) 




2tit Ji(gZ) 

17T/4 / 

(<?)>0 



^-2. £ te^Mj (A.55) 



so that 



$(r, y = 0) ~ RjJ^Yl jl '\ Wl ' n) e-^ r ~ R s \F^C ie -^ (A.56) 

n ' 

Let us now calculate the trapped surface shape in the RS background with 

S(r,y) = RJe~ 2 ^ 1 f°° dqJ (gr) ^f^^ (A.57) 
The continuity condition is 

(V$) 2 + e w (d y $) 2 = (A.58) 
(note the different sign in the exponent, due to the 2-brane RS background). We have 

d y $ = -R s le-^\l l / dqqJo(qr) ^ „ . Q) (A.59) 
Jo hilq) 

(where we have used zl' v {z) = zl v _\{z) — vl v (z)) which means that <9 y $| y= o = 0! And then 

Sg*ly=o = ^ f \{lqfJo{qr)rw{ = R ° l C \<?U<F)Wl (A-60) 

(where we have used zl' v {z) = zl v -\{z) — vl v {z) as well as !§{x) = 2I\(x)/x + 

Again, as before, the integral is zero in perturbation theory, so we must use the contour 
integral as before, and obtain 

v v lm(q)>0 x ' 



We get 



fl2<l>| _ R » l 2 *Vi V .-l/2 3l,nHjl,n) q -„ r 

U y^\y=0 — r- jo 7_jJ\,n / e 

^ V27IT r ^ <fll,n 



R, /2tt/~ _ 2 



-^Vv^" 9ir = - $ l^i (A.62) 
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(since C ± = C 1 (j 1>1 ) i ). 

We are now ready to apply the formalism for the trapped surface in curved background. 
In 2 

^ = f + ay + g V - + ... = C (A.63) 

we have 

a = Ci, g = d 2 y ®\y=o-j(i, f = $\ y =o (A.64) 
We first try a = = (so that ^ = $). Then we need to match 

4 = f + I/W + 9 2 ) + ... (A.65) 
where g = —fqf and /' = —fq~i- We obtain 

C 2 = f(a-y 2 f 1 ) + ... 

4 = / 2 g 2 (l + 0) + ... (A.66) 

which we see that don't match. So we need to put a nonzero a = £i- 
Then we would need to match 

C = f + ay + ... 

4 = /' 2 + a 2 + j/(2a7' + 2 j + 2a<?) + ... (A.67) 

where /' = — q±f and <? = — qff — Aa/l. 

We will try first a = af //, which gives a term in the continuity equation comparable 
with the leading term. We get 

C 2 = / 2 (l + 2^,) + ... 

/ 2 [i+ . 2 r^7» (-^)] + - ( A - 68 ) 



qf + a 2 /P J 1 qf + a 2 /P K I 3 
and by matching the two equations we get 

a = zi±vfEjgg (A . 69) 

however, since j^i = q±l ~ 3.83, there is no real solution! 

Next, we try an order oil/r down from the previous try, namely a = f3f jr. We get 

C 2 = / 2 (l + 2^) + ... 
r 

4 = / 2 g 2 (l + o(^)) + ... (A.70) 
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so now even the r dependence doesn't match. We can see that by adding powers of l n /r n we 
generate again a mismatch of r dependence. 

Conversely, we can try to have an a = (i which is more important than the leading term, 
namely a = afr/l 2 . Then we get 



4 = / 2 ^Al-^) + ... (A.71) 



so again a mismatch of r dependence. 
For higher powers of r n /l n 



f r r 

a = a^( 1 ) n (A.72) 



we need to match 



^ = / 2 [l + 2«4(y) n ] 
r z I 

4 = ^(1-^)1 (A.73) 

and we see that we get a solution for n = 1, a — —3, thus 

f r 2 

a = -SL- (A.74) 

It would seem that we have exhausted all the possibilities, but this is actually not so. 
We can still try 

fir 

a = Ci = a—f (A.75) 
r 



which gives 



VP /3r 

C 2 = f 2 (l + 2a^) + ... 

r 

4 = fqf[l - 2a^ r (/?g 1 + 3a^-)} + ... (A.76) 



rqi Ir 



we can easily see that (3 = —q± makes the two equations equal, so 



-qir 



a = Ci = a / (A.77) 

r 

and strangely a is arbitrary (but probably it will be fixed in a higher order in y). 

In any case, we see that this solution for a does not change the leading order equation at 
y=0, which is 

m * - • /^C ie -*' = 1 (A.78) 
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which has a solution that is approximately 

r H ~ — ln{2q 1 A) 
2q 1 

A = (M)2 27r / C ' 1 2 (A.79) 

The fact that there is a solution can be easily seen by considering the function 

g( r ) = r- Ae~ 2qir => g'(r) = 1 + 2q l Ae~ 2qir > (A.80) 

and the solution we are looking for is given by g(r) = 0. Since g(0) = —A and the function 
is monotonically increasing, it will have a solution. 
At nonzero impact parameter, we get the equation 



■ 2 M 2r 2 * 



* 2 (i-7br) = i ( A - 81 : 

Its solution is thus the zero of the function 



9(r) = r - (1 - ^)Ae- aftr (A.82) 
but now the analysis of the solution is a bit more involved. Indeed, now 

g'(r) = 1 + 2M(1 - ^)e" 2 ^ - ^AT 2 *' (A.83) 

Luckily, we are only interested in the maximum value b max for which g(r) = has a solution. 

If r 2 < b 2 /2, then #(r) > 0. In particular, g(0) = +oo, g(b/y/2) = b/y/2. To have a 
solution of g(r) = 0, we need a minimum, g'{r\) = 0, with g(r{) < 0, so necessarily r\ > b 2 /2. 

But if 1 — b 2 /(2r 2 ) ~ 1 (6 2 /2 significantly lower than r 2 ), the second term in (jA.83|) is 
larger than the third, so g'(r) > (as gi ~ 1/7 3> 1/r). Thus we need instead 



2r 2 



(A.84) 



so if g(ri) < we get 



At b = b max , g(n) = 0, so 



A 

g (n) ~ n - - — e" 25iri (A.85) 



gin 



r 1 = — e"^, 6 2 ~2r 2 



- v / 2-^Zn(4g 1 A) = ^l n [R s q 1 {VK'iCi\^)V2} (A.86) 
2gi gi V ^ 

where R s = y/sG 4 , G 4 = l/(ZMj, 5 ). 
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However, the same physical argument we have used in the AdS case applies. The two 
solutions we have obtained correspond to the trapped surface that would be there if we had a 
four dimensional scattering, and only the function $ would be different (for a = ae~^ ir f /r), 
and the one due to the very large warping outside the 4d IR brane. As the large warping will 
increase the size of the black hole being formed, we have to take the larger trapped surface, 
described by 

f r 2 

a = -3^— (A.87) 
r 

which gives the continuity condition at y=0 



v 2 ' q x l 2 2 V r qj, I 

which has a solution that is approximately 

1 / ( A \ 
ru — ml 

2qi 2qi 



A = { R^ )2 2n ( ^ )2 (A89) 



We can again see that there is a solution by considering the function 

1 



9{r) 



Ar 2 e~ 2qir ~g'{r) = l + {q 1 - -)2Ar 2 e~ 2qir (A.90) 



and the solution for th is given by g(r) = 0. We can see that g'(r) >0ifr>l/gi and 

,1, 1. R 2 18nC 2 , IK . 

g(-) = -[l--±—±) < A.91 

qi qi I 2 4e 2 ji :1 
for sufficiently large R s , so there will be a solution. 
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